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First detection of a VHE gamma-ray spectral maximum from a 
Cosmic source: H.E.S.S. discovery of the VelaX nebula. 
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Abstract. The Vela supernova remnant (SNR) is a complex region containing a number of sources of non-thermal radiation. The 
inner section of this SNR, within 2 degrees of the pulsar PSR B0833-45, has been observed by the H.E.S.S. y-ray atmospheric 
Cherenkov detector in 2004 and 2005. A strong signal is seen from an extended region to the south of the pulsar, within an 
integration region of radius 0.8° around the position (a = 08''35"'00 ! , S = -45° 36' J2000.0). The excess coincides with a 
region of hard X-ray emission seen by the ROSAT and ASCA satellites. The observed energy spectrum of the source between 
550 GeV and 65 TeV is well fit by a power law function with photon index T = 1.45 ± 0.09 s t a t ± 0.2 S y S and an exponential 
cutoff at an energy of 13.8 ±2. 3 s tat ± 4. l sys TeV. The integral flux above 1 TeV is (1.28 ± 0.17 stat ±0.38 sys )xl0"" 11 cnrV. 
This result is the first clear measurement of a peak in the spectral energy distribution from a VHE y-ray source, likely related to 
inverse Compton emission. A fit of an Inverse Compton model to the H.E.S.S. spectral energy distribution gives a total energy 
in non-thermal electrons of ~2 x 10 45 erg between 5 TeV and 100 TeV, assuming a distance of 290 parsec to the pulsar. The 
best fit electron power law index is 2.0, with a spectral break at 67 TeV. 



1. Introduction 



2. Observations with H.E.S.S. 



The region surrounding the Vela supernova remnant (SNR) 
is very well studied across the electromagnetic spectrum and 
contains a number of complex objects, including the SNR 
RX J0852. 0-4622, which has been previously detected in the 
very high energy y-ra y (VHE) range JKatagiri et alJ l2005t 
lAharonian et alJ2005bl). The Vela SNR itself, at a distance esti - 
mated to be ~290 pc (IDodson et alJ2003UCaraveo et al. 2001 ), 
extends over a diameter of ~8°. It is the nearest SNR to contain 
a young active pulsar, PSR B0833-45, with a period of 89 ms 
and a period derivative P of 1.25 x 10~ 13 s s _1 . This implies 
a spin-down luminosity of 7 X 10 36 erg s" 1 and age o f 11,000 
years. Observations by Chandra JHelfand et alJl20o"ll) clearly 
show the torus-like morphology of the compact X-ray nebula 
surrounding the pulsar and have allowed its rotation axis to be 
inferred. The plane of the torus (i.e. pulsar equator) appears to 
intersect the plane of the sky at a position angle of 40.6° relative 
to No rth &N e & Ro manill2004h . and its X-ray spectral index of 
~ 1 .5 ( Mangan oet^Tl2005l) is typical of pre-cooled pulsar wind 
torii. 

The SNR also contains a number of regions of non-th ermal 
emission, including those labelled by Rishbetj] dl958h from 
radio observations as Vela X, Vela Y and Vela Z (which is 
part of the shell of RX J0852.0-46 22). A diffuse emission 
featur e has been detected by ROSAT (Markwardt & Ogelman 
1 19951) in hard X-rays (0.9-2.0 keV), coinciding with the cen- 
tre of the Vela X region. It was first suggested that this fea- 
ture, which is aligned closely with a filament detected at radio 
wavelengt hs, corresponds to the outflow jet from the pole of 
the pulsar llFrail et al.lll997t) . However, the Chandra observa- 
tions showed that this feature lies along the extension of the 
pulsar equator, although bending to the southwest. 

The study of Vela X is important from two perspectives: 
It was the first middle-aged pulsar wind nebula (PWN) to be 
detected, thereby introducing the concept of PWN evolution 
JWeiler & Panagialll980h . Also, it served as the first prototype 
for offset PWN as a result of SNR exp ansion into an inho- 
mogeneous medium (B londin et alj20 0l]). of which G18. 0-0.7 
served as a second example (Gaensl er et al J 120031) . This sec- 
ond PWN was also recen tly identified with VHE source HESS 
J1825-137 JAharonian et alJ2005d) . 

Here we discuss observations of the Vela pulsar and the 
centre of the Vela SNR with H.E.S.S., including the Vela X 
feature. A region of extended y-ray emission is detected, 
to the south of the pulsar position (previously reported by 
iKhelifi et alJ J2005I) . This region was previously o bserved by 
the CANGAROO collaboration jYoshikoshi et alJll997l) and 
emission of VHE y-rays was claimed from a source 0.13° to the 
southeast of the pulsar. Ho wever doubt was late r cast on the sig- 
nificance of this detection JPazelev et alJ20o"i|). A recent p aper 
from the CANGAROO collaboration ( Eno moto et alJ2005l) re- 
tracts the claimed detection, and suggests some evidence for an 
excess coincident with the H.E.S.S. signal. 
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H.E.S.S. is a system of four atmospheric Cherenkov telescopes 
designed to look for astrophysical y-ray e mission above ~100 
GeV; a review of the detector is given bv lHintonl J2004 . The 
Vela region was observed between January and March 2004 
with the complete H.E.S.S. array. In total 10.3 hours of data 
were obtained, after standard data quality selection, at zenith 
angles between 20° and 40°. These observations were taken us- 
ing a method (wobble mode) whereby the source is offset by a 
small angular distance from the centre of the field of view, alter- 
nating between 28 minute runs in the positive and negative dec- 
lination (or right ascension) directions. The observations were 
made at offsets of 0.5° in declination from the position of the 
Vela pulsar, (a = 08 /, 35 m 20 I , 5 = -45° 10' 35" J2000.0), 
which is referred to as position I for the purposes of this paper. 

Following a detection of extended y-ray emission to the 
south of the pulsar in the 2004 dataset, further observations 
were made in 2005 surrounding the position a = 08''35 m 00\ 
6 = -45° 36' (J2000.0), here referred to as position II, which 
was measured as the centre of gravity of the excess. These were 
also taken in wobble mode, and the mean offset for the 2005 
data is 0.5°. After data quality selection 7.2 hours of observa- 
tions were available for analysis, giving total observation live 
time over the two years of 16.4 hours, after dead time correc- 
tion (approximately 10%). The mean zenith angle for the com- 
plete set of observations is 30.2°. 

3. Analysis 

The observations of the Vela region were calibrated using 
standard H.E.S.S. cali bration procedures, as discussed by 
lAharonian et alJ (120041) a nd analysed using the sc aled width 
method, as discussed by lAharonian et al.1 J2005al) . The hard 
cuts, as defined in the latter paper were applied to this source. 
These cuts improve the angular resolution and reduce system- 
atic effects due to uncertainties in the background estimation, 
relative to the standard cuts, albeit at the expense of a higher 
than usual energy threshold (450 GeV for the analysis pre- 
sented here). 

The background level in this analysis was estimated using 
the on-off method, based on a sample of runs taken at simi- 
lar zenith angles (external to the Vela SNR) which contain no 
excess y-ray signal. This method allows the background level 
to be estimated in cases where the integration region for the 
source is comparable to the size of the field of view. The off 
runs are carefully chosen to match well the on source observa- 
tions in zenith angle and system configuration. 

Figure^shows a sky map of the excess in the region south 
of the Vela pulsar. The map has been smoothed with a Gaussian 
function with a sigma of 0.09° and is not corrected for varia- 
tions in the y-ray radial acceptance of the detector, which are 
less than 10% across the integration region. Position II, as de- 
fined above, is marked by a black cross. The extent and posi- 
tion of the excess has been determined by fitting a two dimen- 
sional, elongated Gaussian, convolved with the point spread 
function of the instrument for these measurements, which is 
0.08° (68% containment radius), to an uncorrected excess map 
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Fig. 1. Gaussian smoothed sky map of region surrounding Vela 
pulsar, showing significant emission to the south of the pul- 
sar position, coincident with an X-ray feature seen by ROSAT 
(white contours). The smoothing width used is 0.09°. The con- 
tours corresponding to the strong emission close to the pul- 
sar (Position I) are truncated. The image inset in the bottom 
left corner indicates the size of a point source as seen by 
H.E.S.S., for an equivalent analysis. The solid circle represents 
the H.E.S.S. integration region for the spectral measurement, 
while the dashed circle represents the field of view for the 
ROSAT observations. Position II is marked by a black cross. 
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Fig. 2. Profile of the excess from the extended feature coinci- 
dent with Vela X along minor axis (top) and along major axis 
(bottom); events within 0.8° of the axis are integrated into the 
profile in each case. The two dimensional fit of an elongated 
Gaussian convolved with the H.E.S.S. PSF is also shown as 
a profile for comparison. The position of the Vela pulsar is 
marked by a vertical dotted line in each plot for reference, while 
the integration region in each case for the spectral measurement 
is marked by dashed lines (within 0.8° of position II). 



of the source. The best fit intrinsic width along the major axis is 
0.48° + 0.03°, while the best fit intrinsic width along the minor 
axis is 0.36° + 0.03°. The major axis of the fitted distribution 
is at a position angle of 41° + 7°. The best-fit centre of grav- 
ity of the emission region (a = 08 A 35T, 5 = -45° 34' 40" 
J2000.0) is consistent with position II, within the statistical er- 
ror of 2 arcminutes. The VHE source is thus identified as HESS 
J0835^55. 

Profiles of the excess parallel to, and perpendicular to, the 
major axis of the fit are shown in Figure |2] Profiles of the 
elongated Gaussian fit are shown for comparison. It can be 
seen that, given the statistical errors in each bin, the elongated 
Gaussian fits the excess well, with a^ 2 of 277.8 for 250 degrees 
of freedom. A circular Gaussian fit to the same excess gives a 
width of 0.43° + 0.02°, with a x 2 /dof of 290.0/252. 

The radius of the integration region used for estimating the 
flux and spectrum of the VHE emission was 0.8°, surround- 
ing position II, as defined above. This size was chosen so as 
to encompass the entire excess, as seen in Figure^ A total of 
9610 on-source counts and 7879 o/f -source counts were de- 
tected above 450 GeV, with a normalisation factor of 0.94 be- 
tween the on and off observations, giving an excess of 2152 
events within the integration region. The statistical excess sig- 
nificance is 16. 7c, calculated using the method described by 
iLi & MaHl983h . 



Figure shows the energy spectrum of the signal seen 
from the integration region. The spectrum is well fit (with a 
^ 2 /dof of 13.1/10) by a power law function with photon in- 
dex F = 1.45 + 0.09 s tat ± 0.2sys and an exponential cutoff 
at an energy of 13.8 + 2.3 s tat ± 4.1sys TeV The data may 
alternatively be described using a broken power law function, 
with a spectral break at 13.4 + 3. lstat ± 4.1sys TeV and pho- 
ton index of T\ — 1.7 + 0.1 s tat ± 0.2sys at lower energies 
hardening to Ti = 3.4 + 0.7 s tat + 0.2sys above the break 
energy. The ^ 2 /dof of the broken power law fit, which in- 
cludes a term for the width of the transition region, is 9.0/8. 
A straight power law fit (^- 2 /dof of 74.1/11) is rejected. The 
integral flux above 1 TeV (from the exponential cutoff fit) is 
(1.28+0.17 stat ±0.38 S ys)xl0 -11 crrrV'.The systematic error 
on the flux estimation is estimated at 30%, mainly due to un- 
certainties in the transparency of the atmosphere to Cherenkov 
light. In fitting the energy spectrum for this source, the energy 
of each event is corrected for the detector optical efficiency, 
relative to that used in Monte Carlo simulations to estimate 
the effective area of the instrument. The optical efficiency is 
estimated fro m single muon events de tected during each ob- 
servation run (Lerov et al. 2003; Bolz 2004). The mean energy 
correction is ~20%. 
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Fig. 3. Energy spectrum of y-ray emission from the Vela X re- 
gion. The solid line denotes the best fit of a power law with an 
exponential cutoff. The dashed line represents the best fit bro- 
ken power law spectrum. The bottom panel shows the residuals 
to the exponential cutoff fit. 

In order to test for a point source component in VHE y-rays 
from position I (the pulsar position, as marked in Figure[2j, the 
residual excess within 0. 1° from the pulsar, after subtraction of 
the Gaussian fit to the extended excess, has been measured. No 
significant excess after subtraction is seen, and a 99.9% upper 
limit on the integral flux above 1 TeV, assuming a point source 
at the position of the pulsar, is 7.6 x 10" 13 crrr 2 s _1 . 

4. Discussion 

The new VHE source reported here, HESS J0835-455, is sit- 
uated to the south of the pulsar and the compact X-ray neb- 
ula (as seen by Chandra). The integral flux is estimated to be 
~50% of that of the Crab nebula above 1 TeV. As the distance 
to the pulsar is well measured, one can estimate the size of 
the emission region seen by H.E.S.S. to be 5.1 parsec (full 
length for 68% containment) along the major axis by 3.8 par- 
sec (full width). The luminosity of the emission region in the 
energy range from 550 GeV to 65 TeV can be estimated to be 
L = 9.9 x 10 32 erg s _1 , using the power law fit to the spectrum 
with the exponential cutoff. 

HESS J0835-455 appears to be spatially coincident with 
the X-ray (0.4-2.4 keV) emission as seen by ROSAT (shown 
in Figure 0. It has been suggested jBlondin et al.ll200ll) that 
the Vela X feature corresponds to the pulsar wind nebula, dis- 
placed to the south by the unequal pressure of the reverse shock 
from the SNR. This hypothesis is consistent with the H.E.S.S. 
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Fig. 4. Spectral energy di stribution for the H.E.S.S. and A SCA 
spectral measurements (Markwardt & Qgelmanl Il997l) . The 
two alternative X-ray spectra are described in the text. The fit- 
ted inverse Compton emission (solid line) from the Vela X re- 
gion is shown, given the electron energy distribution described 
in the text. The predicted synchrotron flux is shown for three 
possible magnetic field levels, between 2 fiG and 8 /iG. 

observations which demonstrate conclusively for the first time 
that this feature emits non-thermal radiation. A similar expla- 
nation has been proposed for the emission seen by H.E.S.S. 
(HESS J1825-137) clos e to the pulsar wind nebula Gl 8.0-0.7 
( Aharonian et al . 20053). 

Th e size of the feature described by M arkwardt & Qgelmanl 
d 1995b has been measured as 45 x 12 arcmin 2 , whereas 
the VHE intrinsic size is 58 x 43 arcmin 2 . An extension 
of the X-ray fea t ure to the southwest was suggested by 
iLu & Aschenbachlll2000l) . however this may be unrelated ther- 
mal emission, a suggestion which is supported by the fact that 
no excess is seen in this region by H.E.S.S. 

The larger VHE size probably indicates that synchrotron 
cooling in the X-ray domain was (or still is) important. This 
is confirmed by the detection of a coole d X-ray photon in- 
dex of ~2.1 (Mar kwardt & Qgelmanl 19971 in the emission re- 
gion, which is consistent with the uncooled photon index of 1 .5 
closer to the source. Figure0]shows a spectral energy distribu- 
tion (SED) for Vela X, including the H.E.S.S. result reported 
here. Two possible X-ray spectra are shown, the lower includ- 
ing only the 'Cocoon' feature defined in the above paper, while 
the upper flux includes the emission from the 'SNR' compo- 
nent extrapolated over the H.E.S.S. integration region. 

A simple one-zone model of inverse Compton emission 
from interactions with the cosmic microwave background 
(CMBR) is also shown. The H.E.S.S. observations show the 
expected peak in the IC emission; this result is the first clear 
measurement of such a peak at VHE energies. The total energy 
content in non-thermal electrons between 5 TeV and 100 TeV is 
2.2 x 10 45 erg, for a distance of 290 parsec. The best fit electron 
power law index is 2.0, with a break at 67 TeV and a post-break 
index of 9. The ^ 2 /dof of this fit is 10.2/9. In Figure |4] three 
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possible synchrotron spectral energy distributions are shown, 
for magnetic field levels in the emission region from 2 fiG to 
8 fiG. The synchrotron distributions are not constrained by the 
X-ray data, a combined fit on the X-ray and y-ray data would 
more closely reproduce the measured X-ray spectral slope. A 
magnetic field of this magnitude may be expected in the sce- 
nario of a nebula displaced from the pulsar by an asymmetric 
shock. 

This result demonstrates the usefulness of H.E.S.S. obser- 
vations in clarifying the complex morphology of this source. 
VHE observations of inverse Compton scattering of the CMBR 
allow direct inference of the spatial and spectral distribution of 
non-thermal electrons in a PWN, independent of contamination 
by thermal emission or variations in the local magnetic field. 
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